Towards testing interacting cosmology by distant type la 

supernovae 

Marek SzydlowskiJ*] Tomasz Stachowiak, and Radoslaw Wojtak 
Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland 

(Dated: February 3, 2008) 

Abstract 

We investigate the possibility of testing cosmological models with interaction between matter 
and energy sector. We assume the standard FRW model while the so called energy conservation 
condition is interpreted locally in terms of energy transfer. We analyze two forms of dark energy 
sectors: the cosmological constant and phantom field. We find a simple exact solution of the models 
in which energy transfer is described by a Cardassian like term in the relation of H 2 (z), where H is 
Hubble's function and z is redshift. The considered models have two additional parameters (fimt) n ) 
(apart the parameters of the ACDM model) which can be tested using SNIa data. In the estimation 
of the model parameters Riess et al.'s sample is used. We also confront the quality of statistical fits 
for both the ACDM model and the interacting models with the help of the Akaike and Bayesian 
informative criteria. Our conclusion from standard best fit method is that the interacting models 
explains the acceleration of the Universe better but they give rise to a universe with high matter 
density. However, using the tools of information criteria we find that the two new parameters 
play an insufficient role in improving the fit to SNIa data and the standard ACDM model is still 
preferred. We conclude that high precision detection of high redshift supernovae could supply data 
capable of justifying adoption of new parameters. 
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I. INTRODUCTION 



The main aim of the paper is to use the distant type la supernovae (SNIa) compiled by 
Riess et al. Q] to analyze the interacting cosmology (or cosmology with energy transfer) 
which was recently proposed by Zimdahl and Pavon [2| and Zhang 3| as a way to explain the 
coincidence conundrum of the ratio of dark matter to dark energy density at present. We 
show that allowed intervals for two independent parameters which characterize interacting 
cosmology give rise to high matter density, which can hardly be reconciled with the current 
value obtained from cluster baryon fraction measurements. 

The major development of modern cosmology, which mainly concentrates on observational 
cosmology, seems to be a discovery of the acceleration of the universe through distant SNIa. 
In principal there are three different theoretical types of explanation of acceleration. In the 
first group the existence of matter of unknown forms violating the strong energy condition 
is postulated (dark energy). In the second group instead of dark energy some modification 
of the FRW equation is proposed. Between these two propositions lies an idea of interacting 
cosmology in which the source of acceleration lies in interaction between dark energy and 
dark matter sector. 

The starting point of any investigations in observational cosmology is the naive FRW 
model (the ACDM model) in which matter and dark energy sectors do not interact with each 
other. The next step in making this model more realistic is to include the interaction between 
dark matter and dark energy. We start with the FRW standard cosmological equations: 
the Raychaundhuri equation and conservation condition. They assume the following forms 
respectively 

- = ~(p + 3p) (i) 

a o 

P = SH(p + p), (2) 

where a is the scale factor, p and p are energy density and pressure of perfect fluid, a source 
of gravity; H = a/a is the Hubble function, and a dot means differentiation with respect to 
the cosmological time t. It would be convenient to rewrite relation (J2J) for mixture of dust 
matter and dark energy 
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where the total pressure and energy density are 



p = + w x px, w x = Px/px = const (4) 
P = Pm + Px (5) 

The standard postulate of the ACDM model is that both terms in eq. (J3J) are equal to zero, 
or that the matter components evolve independently. We, however, assume the existence 
of interaction introducing a non- vanishing function j(t) such that eq. © is still satisfied 
locally 

~M 3 )=7(*) (6) 
1 (Pxa 3 ^) = - 7 (t) (7) 



a 3(l+^x) dt 

Note that 7(£) in ©-© has the interpretation of the rate of change of energy in the unit 
comoving volume. 

We only assume the existence of the energy transfer without specifying its physical mech- 
anism, then test it by astronomical observations. The idea of such decomposition of the 
conservation condition ©-© is not new Q. The interacting quintessence models with 
pressureless cold dark matter fluid can give rise to the scale factor relation px/Pm = cl^ 
P, 0|, where £ is a constant parameter. This relation was also tested by CMB WMAP 
observations and statistically investigated in the context of SNIa observations [?J but with 
the help of a different exact solution. In the paper of one of us , this idea was investigated 
independently without any physical assumption or about the form of such an interaction or 
an exact solution. Under such general assumptions, that the FRW dynamics is valid, the 
existence of energy transfer to the dark energy sector from the dark matter section is favored 
by SNIa data. However, if we apply the prior f2 mj o = 0.3 as claimed by independent galactic 
observations, then the transfer between dark matter and phantom sector is required. The 
proposition of different physical mechanism of e nergy e xchange from dark energy to dark 
matter sectors was proposed by many authors jj], lloillll| . 

We consider models where the rate of energy transfer is proportional to the rate of change 
of the scale factor ^(t) = f(a)H. 

While the standard statistical methods of maximum likelihood used for estimation of the 
model parameters favor the energy transfer between dark matter and phantom dark energy 
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sector, we demonstrate by the means of the Akaike and Bayesian information criteria that 
the ACDM model with interaction is selected. 



II. THE FRW MODEL WITH INTERACTION 



We assume for simplicity of presentation without loss of generality a simple flat model 
for which we have the Friedmann first integral in the form 

3d 2 = a 2 (p m + p x ) (8) 

In the case of non-flat models relation (JSJ) is of the form 

3a 2 = a 2 (p m + p x )-3k (9) 

We add to equation (jSJ) the generalized adiabatic condition 

j t p m = -3H Pm + 1 (t), (10) 

j t p x = -3(1 + w x )Hp x - T (t), (11) 

J t H^-H 2 -\ {Pm + p x ) (12) 



from which we obtain 



!^ s/ . taJ = ^ + JL. (i3) 

a(ma) -np m 

*£H± S I M = -2- (14) 
a{ma) Hp x 

*mpJk)_ 3wx + ul + ±) (15, 



d(lna) H \p m p x J 

where I a (Pm) and I a {px) have the interpretation of elasticity (logarithmic slope) of the energy 
densities with respect to the scale factor. If energy transfer vanishes then the corresponding 
elasticity function becomes constant. From (|T3*j) -()14 |) we can simply derive the relation 



„ \ d In & I O 0/ TT 



\px J d(lna) p m p x 



Therefore scaling solutions Pr m/px_p^ cl^ can only be realized if a very special form of param- 

110,13], namely 

7 (t) = lQ H^-^ = lop-'H, (17) 



eterization ^(t) is assumed |4|,|y, l2(, namely 



Pm PX 



where p is the harmonic average of p m and px- It is interesting to discuss the existence 
of scaling solutions for other parameterizations of 7(£). For this purpose the methods of 
dynamical systems can be useful because they offer possibility of investigating all admissible 
evolutional paths for all initial conditions. Let us consider the following distinguished case. 

Let 7(t) = 7o/(a)if and f(a) = a m . In this case one can find exact solutions determining 
energy densities 

p m = ^- m+2 + ^ (18) 
m + 5 or 

„ 70 m+2 , Aa,0 /, q x 

P* m + 5 + 3wx a + a 3(l+^) ^ 

where in the special cases of m = —5 and m = — 5 — 3wx instead of a power law function 
a logarithmic function appears; A m and Ax,o are constants. Note that the contribution to 
the total energy arising from interaction is of the type a~ m ~ 2 , therefore it can be modelled 
by some additional fictitious fluid satisfying the equation of state p- mt = — (m + 5)pi nt /3. 
In terms of redshift z we obtain the basic relation 

H 2 = Hi (Q int (l + z)~ m - 2 + fi m , (l + zf + fi x ,o(l + zf 1+w ^ + Q M (1 + z) 2 ) , (20) 

where 

Qint,o = 7 — X r7o, (21) 

[m + 5)(m + 5 + 3wjs:) 

and the index zero means that quantities are evaluated at the present moment. The relation 
will be useful in the next section where both Q- int and m parameters will be fitted to SNIa 
data. 

For formulation the corresponding dynamical system it would be useful to represent this 
case in the form of the Hamiltonian dynamical system. For this case it is sufficient to 
construct the potential function 

V = (22) 
6 V ; 



which build the Hamiltonian 



a 2 



H = - + V(a). (23) 



We find that the Hamiltonian system is determined on the zero level Ji — and the potential 
function in terms of the scale factor expressed in units of its present value ao, x — is 



V{x) = - l -tt mfl x- 1 - ]f wxfi x-^ 1+w ^ +2 - l -n- int x~ m - ±Q kt0 . (24) 



The Hamiltonian constraint gives rise to the relation Yli = 1- The corresponding dy- 
namical system is of the form 



x = y (25) 

»--£ < 26 ' 

and Ti, = plays the role of first integral of ([25 )1 - ([26 )1 . Advantages of representing dynamics 
in the form ()25j1 - ()2f)j) is that the critical points as well as their character can be determined 
from the geometry of the potential function only. Only the static critical points, admissible 
in the finite domain, y = 0| zo , ^ = are saddles or centers. The eigenvalues of the 
linearization matrix are 

a - = ± fWT n <27) 

The phase portraits on the phase plane (x,x) are shown in Fig. [I] (w — —1, "decaying" 
vacuum cosmology) and in Fig. |21 (w = — | , phantom fields) . We can observe the existence 
of a single static critical point located on the x-axis and the topological equivalence of these 
phase portraits in the finite domain. In both cases we have static critical points situated 
on the x-axis. They are of saddle type because the potential function is upper convex. We 
have also the critical points at infinity which represent stable nodes. They represent the de 
Sitter model in Fig. ^ and the big rip singularity in Fig- El The big rip singularity is generic 
future of the model violating weak energy condition at which the scale factor is infinite at 
some finite time [l^ . 

In connection with the present value of the density parameters f2 mj o and flx,o there arises 
the coincidence problem. Why are the energy densities of dark energy and dark matter of 



the same order amplitude at the present epoch? 
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14, 



15. 161. In this context we can 



use some interesting idea formulated by Amendola and Tocchini-Valentini [17[ who argue 
that the present state of the accelerating universe is described by the global attractor in 
the phase space. It seems to be reasonable to treat the evolution of the interacting universe 
model with tools of dynamical system methods. From the physical point of view the critical 
points represent asymptotic state of the system reached in infinite time. Let us assume that 
7(t) is proportional to the Hubble parameter, i.e. cx j(p m , Px)H WWW then after 
the reparametrization of the time variable t t: Hdt = dr, we obtain the two-dimensional 
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FIG. 1: The phase portrait for the A interacting cosmological model. The trajectories of flat model 
separate the domains of closed and open models 



closed dynamical system in the form 

dfm _ 
dr 

dpx 
dr 



-3p m + 7(p m; PX) 
-3(1 + w) -7(Pm,Px)- 



(28) 
(29) 



System (j2~8*j) - (}2T)j) has a critical point at which 



^x,o/^m,o = -(1 + w) 1 - 

In the case of phantom (w = —4/3) we obtain Qx,o/^m,o — 3. This critical point is stable if 
a trace of a linearization matrix is negative at this point, i.e., <9(— 7)/<9p m > d(—^f)/dpx — 2. 
If 7 oc —p m px it means that Qx,o > ^m,o- 

In Fig. El we observe that the most probable values of n — —m — 2 and f2i n t,o are close to 
zero. In Fig. 0] there is two disjoint region which are preferred, while values of n and r2 int 
close to zero are excluded at the 95% confidence level. 



future big rip- 
singularity 



Einstein 
static solution 




FIG. 2: The phase portrait for the phantom interacting cosmological model. 

For deeper statistical analysis it is useful to consider the probability distribution function 
(PDF) over model parameters. On Fig. |SJ |HJ El the one-dimensional PDF are demonstrated 

for fimt.O, n > ^A,0- 

III. DISCUSSION OF RESULTS 



We use the maximum likelihood method to estimate two basic parameters of the cosmol- 
ogy with interaction. The priors we assume for numerical computations are: — 10 < n < 10, 
— 1 < Qi n t < 1, < fi m) o < 1 (or fixed at 0.3), and the other density parameter determined 
from £\ Qi — 1. From this statistical analysis we obtained that with assumed fl m ,o = 0.3 the 
interaction between "decaying" A and matter sectors is negligible while the interaction be- 
tween phantom and matter sector is required (see Table |1J). The analogous analysis without 
any priors contains (see Table [H} . 

To select between these models we use the Akaike and Bayesian information criteria 
(AIC and BIC) Q 19 ,Q These criteria were proposed to answer the question whether 
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FIG. 3: Confidence levels on the plane (f2; nt ,n) for the interaction between matter and the A 
vacuum energy sectors; the prior on O mi o = 0.3 is assumed. 

or not new parameters are relevant in the estimation. In the generic case introducing new 
parameters allows to improve the fit to the dataset, however information criteria try to 
measure whether the improvement in fit is large enough to consider it statistically significant. 

When we consider the interacting cosmological model we add two new parameters de- 
scribing the interaction to the base ACDM model. Using the AIC and BIC we compare the 
fits of these two models to the most recent sample of SNIa (see Table ITTij) . 

To improve the statistical analysis and obtain more precise physical results it is necessary 
to extend the sample to higher redshifts z > 1.5. Given two mechanisms of interaction 
between dark energy and dark matter sectors, we found that the ACDM model is favored 
over the "decaying" A interacting model (AINT) but it is disfavored over phantom interacting 
model (PhINT) for both criteria. The same result is obtained independently of using the 
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FIG. 4: Confidence levels on the plane (f2i n t> n ) f° r the interaction between matter and the phantom 
sectors; the prior on Q m q = 0.3 is assumed. 
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FIG. 5: The PDF for (fiint,o) when matter interacts with the "decaying" vacuum (left) and matter 
interacts with phantom field (right); the prior on O m) o = 0.3 is assumed. 
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FIG. 6: The PDF for n when matter interacts with the "decaying" vacuum (left) and matter 
interacts with phantom field (right); the prior on ^ m ,o = 0.3 is assumed. 
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FIG. 7: The PDF for (Oa,o) when matter interacts with the "decaying" vacuum (left) and matter 
interacts with phantom field (right); the prior on £l nij o = 0.3 is assumed 

prior on fi m ,o- 

It is an interesting question whether it is possible to determine the direction of energy 
transfer from SNIa data. To this aim we calculate the probability that 70 is greater than 
zero -P(7o > 0) as well as find the most probable values of the parameters. We have 



(30) 



P( 7o > 0) = 35% for AInt models 

P( 7o > 0) = 82.81% for Phlnt models, 

under the constraint f2 m = 0.3. Therefore, with the confidence level of 95% we can't exclude 
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TABLE I: Results of the statistical analysis of the model obtained both for Riess et al. samples 
from the best fit with minimum \ 2 (denoted with BF) and from the likelihood method (denoted 
with L). The same analysis was repeated with fixed £l m ,o = 0.3 (denoted F). 



sample 


" "m,0 


S l\r,i- n 71, — 
u inx,u ' v 


— 171 — 2 


^A,0 M 


x 2 


method 


ACDM 


0.31 






0.69 15.955 


175.9 


BF 




0.44 






0.69 15.955 




L 




F 0.30 






F0.70 15.955 


175.9 


BF 




F 0.30 






F0.70 15.955 




L 


AINT 


0.43 


0.28 


-10.00 


0.29 15.905 172.1 


BF 




0.44 


0.34 


-2.70 


0.20 15.935 




L 




F 0.30 


0.07 


-10.00 


0.63 15.925 


175.2 


BF 




F 0.30 


0.00 


0.00 


0.70 15.945 




L 


PhINT 


0.46 


0.25 


-10.00 


0.29 15.905 


172.2 


BF 




0.44 


0.26 


-2.60 


0.62 15.935 




L 




F 0.30 


0.03 


5.10 


0.67 15.935 


173.5 


BF 




F 0.30 


0.03 


0.00 


0.67 15.945 




L 



any direction of energy transfer. Because this result is inconclusive, we check the sign of 70 
for the most probable values of the parameters. In both cases, we obtain (independently 
of assuming the fixed prior for fi m ,o) that 70 < 0. Let us also note that in the case of the 
PhlntCDM model the values fii n t,o = and n = are excluded with the confidence level of 
95%. 

Both the maximum likelihood method and the best fit method are used principally to 
estimate the values of parameters. Additionally the best fit method allows to choose the 
best model when a few models are to consider. However there is a serious shortcoming of the 
best fit method because models with greater number of parameters are generally preferred. 
To overcome this drawback and push the statistical analysis one step further we employ the 
model selection criteria. It can be seen as the second stage of statistical analysis after the 
estimation. Both of the used information criteria, the AIC and BIC, put some penalty on 
adding additional parameters assigning a threshold which should be exceeded to consider a 
model with a greater number of parameters to be a better fit. 
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TABLE II: Model parameter values obtained from the minimization procedure carried out with 
the Riess et al. sample. 



sample 




^int,0 


^A,0 


n = —m — 2 


ACDM 


0.31+g-gt 

— U.U4 




0.69+^4 

— U.U4 






0.30 




0.70 




AINT 


0.30 


n S4+ ' 39 

^■^-0.33 

0.001°;?? 


70+ - 17 


n 7n +i.80 

o oo +1 - 50 


PhINT 


n 44+ - 10 
0.30 


26+ ' 50 
U - ZO -0.46 

01+ - 43 


R9+ - 24 
u - DZ -0.76 

67+ - 03 
U - D/ -0.43 


Z - DU -3.10 

o oo +1 - 50 



TABLE III: The values of AIC and BIC for distinguished models both for ri m o fitted and for 
^m,o = 0.3. 



case 


AIC (fitted Q mfi 


) AIC (fi mi0 = 


0.3) BIC (fitted O m , ) 


BIC (n m>0 = 0.3) 


ACDM 


179.9 


177.9 


186.0 


181.0 


AINT 


180.1 


181.2 


192.3 


190.4 


PhINT 


180.2 


179.5 


192.4 


188.7 
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